Air Entrainment by Viscous Contact Lines 
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The entrainment of air by advancing contact lines is studied by plunging a solid plate into a 
very viscous liquid. Above a threshold velocity, we observe the formation of an extended air film, 
typically 10 microns thick, which subsequently decays into air bubbles. Exploring a large range of 
viscous liquids, we find an unexpectedly weak dependence of entrainment speed on liquid viscosity, 
pointing towards a crucial role of the fiow inside the air film. This induces a striking asymmetry 
between wetting and dewetting: while the breakup of the air film strongly resembles the dewetting 
of a liquid film, the wetting speeds are larger by orders of magnitude. 



Objects that impact on a liquid interface can entrain 
small bubbles of air into the liquid. This happens for ex- 
ample when raindrops fall in the ocean [if or when liquid 
is poured into a reservoir at sufficiently large speeds Q ■ 
Such entrainment of air is often a limiting factor in indus- 
trial applications such as coating and nano-scale printing 
techniques, where the bubbles disturb the process 0, 
A well studied case is the entrainment of air by very 
viscous lets impacting on a reservoir of the same liq- 
uid 0, li , 0] • The onset of entrainment is essentially 
determined by the properties of the liquid, Ue ^ 7/??£, 
which reflects a balance of the liquid viscosity rji and the 
surface tension 7. Changing the nature of the gas only 
has a minor influence on the entrainment process Q . 

A very different picture has emerged recently in the 
context of drops impacting on a walL for which the pres- 
ence of air has a dramatic effect |8l4lOl|. It was found 
that splashing can be suppressed completely by reducing 
the air pressure to about a third of the atmospheric pres- 
sure. This caused huge excitement, in particular because 
such a pressure reduction does not lead to any change 
of the gas viscosity rjg 11- 13]: pressure only affects the 
gas density, and thus the speed of sound and the mean 
free path in the gas. A similar paradox is encountered 
for air entrainment by rapidly advancing contact lines, 
where a liquid advances over a surface that it partially 
wets [3, [l^l8l | . Once again, it was found that depressur- 
izing the gas leads to a signiflcant increase of the thresh- 
old of air entrainment 18 , 3] ■ This contradicts the clas- 
sical viewpoint that, for given wettability, the contact 
line speed depends mainly on the liquid properties as 
^ llVi (Hi 13" 22 1: with minor influence of the gaseous 
phase. 

In this Letter we reveal the role of the air for advanc- 
ing contact lines in a paradigmatic system: a partially 
wetting solid plate is plunged into a reservoir of viscous 
liquid (Fig.[T]). At sufficiently high speeds we observe the 
entrainment of an air film, typically 10 microns thickness, 
which subsequently decays into bubbles. The liquid vis- 
cosity is varied over more than two decades by using 
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FIG. 1: Air entrainment by a contact line when a solid plate 
is plunged into a viscous liquid. (a,b) Sketch of dynami- 
cal meniscus below the threshold of air entrainment (a), and 
above the threshold of air entrainment (b). In the latter case 
a film of thickness h is formed. (c,d) Experimental images of 
the entrained air film in silicone oils of varying viscosity r^i. 
(c) For r]i = 0.1 Pa.s, an extended air film is entrained behind 
the contact line moving downwards into the bath. The film 
is destroyed by the formation of "rewetting" holes. The di- 
agonal dark line is the refiection of a straight wire that gives 
an impression of the interface profile, (d) For a more viscous 
oil, Tji = 5 Pa.s, the entrained air film is rapidly destabilized 
to form small bubbles (dark spots). 



silicon oils of different molecular weights. It is found 
that the entrainment speed Ue changes much less than 
the expected scaling l/?7£- Using an approximate hy- 
drodynamic model we argue that this can be attributed 
to the flow of air into the strongly confined film, mak- 
ing the contact line velocity strongly dependent on both 



gas and liquid viscosities. This induces a striking asym- 
metry between wetting and dewetting: the same hquid 
can advance much faster than it recedes, by orders of 
magnitude. 

The experiments are carried using sihcon oils (PDMS, 
Rhodorsil 47V series) with dynamic viscosities rje = 
0.02,0.10,0.5,1.0 and 5 Pa.s. These liquids are essen- 
tially non-volatile, insensitive to contamination, while 
the surface tension 7 = 22 mN.m~^ and density p = 
980 kg.m"'^ are approximately constant for all viscosities. 
The reservoir is a transparent acrylic container of size is 
29x15x13.5 cm, which is much larger than the capillary 
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1.5 mm. The substrate con- 



sists of a silicon wafer (circular, diameter 10 cm), which 
is coated by a thin layer of fluorinated material (FC 725 
(3M) in ethyl acetate). For all liquids, this results into 
static contact angles between 51° and 57°. The wafer 
is clamped onto a 10 mm thick metallic blade screwed 
to a 50 cm long high-speed linear stage. The combi- 
nation of using controlled speeds and very viscous liq- 
uids one avoids complexities of S2lashiiig as well as the 
formation of interface cusps 0, li, 16-3]. In addition 
the effect of inertia is eliminated both in the gas and in 
the liquid; the relevant Reynolds numbers are at worst 
~ 1, but typically orders of magnitude smaller. For 
each liquid, we plunge the wafer into the reservoir at 
different plate velocities Up, up to 0.7 m/s. The pro- 
cess is recorded using a high-speed Photron SA3 camera 
(2000Hz, 1024x1024 pixels). The plate and contact line 
velocities are then extracted from space-time diagrams 
using a correlation technique with a sub-pixel resolution, 
leading to a precision better than the percent. Repro- 
ductibility is achieved within 15%. The film thickness h 
is determined from an accurate measurement of the vol- 
ume of air entrained in the bath. The air bubbles trapped 
at the end of an experiment are imaged with Nikon D300s 
(4288x2848 pix) mounted with macro-lens. The size of 
each bubble is determined automatically within a 10 /im 
resolution, by fitting an elliptic shape. Most of the un- 
certainty on h results from the estimate of the surface 
covered by the air film. 

The experimental scenario is presented in Fig. [T] At 
small speeds we observe that the contact line equilibrates 
to form a stationary meniscus and no air is entrained 
into the liquid (Fig. [l}i). Above a critical velocity, how- 
ever, the contact line keeps moving downward into the 
reservoir and deposits a thin film of air (Fig. [TJd) . Subse- 
quently, the air film breaks up into small air bubbles. The 
dynamical structures that appear during the breakup of 
the film turn out to depend strongly on the oil viscosity, 
as can be seen from the experimental images of Fig. [TJ;d. 
For the most viscous liquids the air films are very frag- 
ile and rapidly break up into smaller bubbles (Fig. [TJl). 
However, for the lower viscosities one observes the forma- 
tion of very extended air films (Fig. [TJ;). At the front of 
the film, the contact line develops a ridge-like structure 




FIG. 2: "Rewetting" holes (air film invaded by liquid) as 
the inverse of the classical dewetting of a liquid film. (a,b) 
Rewetting: once the silicon oil reestablishes contact with the 
solid, one observes the growth of a circular zone that invades 
the air film. The moving front collects the air in a thick rim 
that is clearly visible in the image. Typical contact line speeds 
range from 1 to 10 cm/s for the silicone oils used in this study. 
(c,d) Dewetting: a thin film of silicon oil dewets in the form 
a circular hole. Note that for the same liquids as in (a) the 
typical dewetting speeds are smaller by orders of magnitude 
(from 40^m/s to Icm/s). 



that is common for dewetting of liquid films |23l - l25| . The 
peculiarity of the present experiment is that in this case 
the air is dewetted, not the liquid. 

An even more striking analogy with classical dewetting 
of liquids is the nucleation of nearly circular regions in- 
side the film (Fig.[T]:). However, the circles now represent 
regions of rewetting, where the liquid reestablishes the 
contact with the solid. These "rewetting holes" can be 
considered as the inverse of the "dewetting holes" , since 
the roles of air and liquid are exchanged. Figure [5^ shows 
a close-up of a rewetting hole (cross-section sketched in 
Fig. [5]d). The radius of the holes increases linearly with 
time, and the advancing contact lines collects the air in- 
side a thick rim. While this is analogous to the inverse 
problem of the dewetting holes (Fig. Htd), the process is 
by no means symmetric: the rewetting holes grow with a 
velocity Ue that is orders of magnitude faster than their 
dewetting counterparts, up to a factor 1000 for the liquids 
used in this study. 

We further quantify the velocity of air entrainment for 
different liquid viscosities. A first measurement of Ue 
is obtained from the growth velocities of the rewetting 
holes as in Fig. [5] A second entrainment velocity can 
be extracted by a selecting a central part of the front of 
the air film, for which we obtain the contact line veloc- 
ity Uci in the frame of the liquid reservoir. Figure [3^ 
reports the measured values for rjf — 0.1 Pa.s, showing 
that the contact line velocity increases linearly with plate 
velocity Up. Interestingly, however, the relative velocity, 
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FIG. 3: (a) Contact line velocity Uc (triangles) as a func- 
tion of the plate velocity Up, for 77^=0. IPa.s. The relative 
velocity Up — Uc (circles) is independent of the plate speed, 
(b) Entrainment speed Ue for different liquid viscosities rji, 
measured in two ways: the circles represent the relative ve- 
locity with respect to the plate, Up — Ud, while the triangles 
are the rewetting speeds of the holes as in Fig. [J] (c) Film 
thickness h for different rje, taken at constant plate velocity 
Up = 0.67 m.s"\ The open circle was taken for a liquid jet 
of glycerol entraining air at the same speed 0]. 



Uci, turns out to be independent of the plate ve- 



locity and therefore seems to be an intrinsic property of 
the advancing contact line. This is why we may con- 
sider the growth velocity of the rewetting holes to be 
an independent measurement of the entrainment veloc- 
ity. The resulting entrainment velocities Ue are shown 
in Fig. ^jp, as a function of liquid viscosity iji. Indeed, 
the two experimental definitions of Ue agree very well for 
the smallest rjg (closed circles are based on the front of 
the film, open triangles correspond to rewetting holes). 
For larger rji, the film rapidly destabilizes and it is more 
difficult to define the front of the film. This induces a dif- 
ference between the two types of velocity measurements 
of about a factor 2; the hole velocities are certainly more 
reproducible in this very viscous regime. 

The key result of the velocity measurements is that, al- 
though Ue decreases with liquid viscosity, the dependence 
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FIG. 4: Dimensionless entrainment speed, Ca = Uerje/'y, ver- 
sus the viscosity ratio rjg/rje for: silicon-oil/air (•), silicon- 
oil/air (■ after IS]) and various liquids/air (A after 26]). 
Oil slip lengths \i = 10~^£-y, and air slip lengths Xg = 10~'^£^ 
(solid line) and 10~^£^ (dotted line), which correspond to 
mean free paths £mf ~ 70 nm (solid) and 7 /im. 



is clearly much weaker than the expected ^ l/?7^- The 
entrainment speed is reduced by a factor 10, while vis- 
cosity is varied by a factor 250 (Fig.tSj)). Since the liquid 
inertia is negligible for these highly viscous liquids, this 
means that the properties of the air must have a signif- 
icant infiuence on the entrainment speed. On the other 
hand, the speed is not determined by the air alone, since 
that would yield no dependence on rji at all. To reveal 
the interplay between air and liquid phases, we introduce 
a dimensionless capillary number, Ca — UeTji/^, that is 
based on the liquid viscosity. The experimental results 
are represented in Fig. 21 showing Ca versus the ratio of 
gas and liquid viscosities rig/rjg (closed circles). Clearly, 
the capillary number for air entrainment displays a de- 
pendence that is much stronger than ~ ln(77£/77o), which 
is the scaling for air entrainment by liquid jets [H, 0] and 
the prediction by Ref. [i^l- The air thus has a much 
larger influence than expected. On the same figure we 
collected data from the coating literature, based on tapes 
running continuously in a bath, showing a similar trend 
(various symbols, see caption). Note that in these exper- 
iments the contact line typically develops a sharp cusp 
from which small air bubbles are emitted, rather than an 
extended air film. 

In order to understand the influence of air, one can 
extend the common lubrication approximation to large 
slopes and 2 phase flow, but small curvature of the in- 
terface. In the spirit of 28, [2^, one can start from 
the analytical solutions of flow in a wedge due to Huh 
& Scriven 30j], as detailed in the Supplementary Mate- 
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FIG. 5: Insets: shape of the dynamical meniscus at the en- 
trainment threshold for rig/rji = 10~*, as predicted by the 
hydrodynamic model [13]. Main graph: local slope S as a 
function of the local thickness h. Lengths are in units of l-y. 

rial {i?!. Numerical solution of the model provides the 
provides the capillary number for entrainment Cag shown 
in Fig. 0] as the solid line. This captures the order of mag- 
nitude for the entrainment speed as well as the unexpect- 
edly strong dependence on viscosity ratio. However, due 
to the strong curvature of the interface, the approxima- 
tion cannot be expected to give a fully quantitative pre- 
diction. Within this model, the viscous dissipation rates. 
Eg and Ei in the gas and the liquid can be estimated as 
a function of the local slope 0: 
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(1) 



which is slightly different from the estimates in [IJ, [15[ . 
Clearly, the dissipation in the gas can dominate over liq- 
uid dissipation for contact angles very close to tt, with a 
cross-over occurring when rig/rii ^ [tt — 0)^. As shown 
in Fig. [51 such angles indeed appear naturally when the 
liquid advances rapidly, while this is not the case when 
the liquid recedes. This explains the marked asymme- 
try between wetting and dewetting: the faster advancing 
speed is due to the reduction of dissipation near 9 k, tt. 

In this Letter we experimentally showed that the en- 
trainment speeds of advancing contact lines do not scale 
as 7/»7^, but has a much weaker variation with liquid 
viscosity. We explain this by the influence of the air 
flow when the local angle of the interface is close to tt. 
Can such a scenario explain the observed increase of en- 
trainment speed when depressurizing the air? A pres- 
sure reduction does not affect the gas viscosity, but it 
does increase the mean free path by a factor ~ Patm/p- 
Since the mean free path under ambient conditions is 
approximately 70 nm, it is pushed well into the micron 
range when pressure is reduced by a factor 100 



m 1 



The mean free path then becomes comparable to the 
film thickness measured experimentally (Fig. [3};). In- 
terpreting the mean free path of the gas as an effective 
slip 271, |31| , this indeed yields an increase of Cae (dotted 
line in Fig. S]). This suggests that depressurized air acts 
as a Knudsen gas when confined close to the contact line. 
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